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ABSTRACT 

Observations from the Hopkins Ultraviolet Telescope and the Hubble Space 
Telescope have recently detected Hell absorption along the lines of sight to two high 
redshift quasars. We use cosmological simulations with gas dynamics to investigate 
Hell absorption in the cold dark matter (CDM) theory of structure formation. We 
consider two 17 = 1 CDM models with different normalizations and one open universe 
(Qq = 0.4) CDM model. The simulations incorporate the photoionizing UV background 
spectrum computed by Haardt & Madau (1996), which is based on the output of 
observed quasars and reprocessing by the Lya forest. The simulated gas distribution, 
combined with the Haardt & Madau spectral shape, accounts for the relative observed 
values of thi and THeii, the effective mean optical depths for HI and Hell absorption. 
If the background intensity is as high as Haardt & Madau predict, then matching 
the absolute observed values of thi and rneii requires a baryon abundance larger (by 
factors between 1.5 and 3 for the various CDM models) than our assumed value of 
Qbh 2 = 0.0125. The simulations reproduce the evolution of THeii over the observed 
redshift range, 2.2 ^ z <^ 3.3, if the Hell photoionization rate remains roughly constant. 

Hell absorption in the CDM simulations is produced by a diffuse, fluctuating, 
intergalactic medium, which also gives rise to the HI Lya forest. Much of the Hell 
opacity arises in underdense regions where the HI optical depth is very low. We 
compute statistical properties of the Hell and HI absorption that can be used to test 
the CDM models and distinguish them from an alternative scenario in which the Hell 
absorption is caused by discrete, compact clouds. The CDM scenario predicts that a 
substantial amount of baryonic material resides in underdense regions at high redshift. 
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Hell absorption is the only sensitive observational probe of such extremely diffuse, 
intergalactic gas, so it can provide a vital test of this fundamental prediction. 

Subject headings: quasars: absorption lines, Galaxies: formation, large-scale structure 
of Universe 

1. Introduction 

The Lyman-a forest seen in quasar spectra (Lynds 1971; Sargent et al. 1980) is produced by 
absorption from diffuse hydrogen at high redshift. The big bang model predicts that approximately 
25% of primordial baryonic matter should also be in the form of helium, and that there should 
therefore be corresponding absorption at shorter wavelengths from the Lya transition of singly 
ionized helium (Hell). It has recently become possible to detect redshifted Hell absorption with 
space-based, ultraviolet (UV) observations (Jakobsen et al. 1994; Davidsen, Kriss, & Zheng 1996; 
Hogan, Anderson, h Rugers 1996). Recent hydrodynamic cosmological simulations of cold dark 
matter (CDM) models have been remarkably successful in reproducing the observed properties 
of the HI Lya forest (Cen et al. 1994; Zhang, Anninos, & Norman 1995; Hernquist et al. 1996, 
hereafter HKWM; Miralda-Escude et al. 1996; Dave et al. 1997; Zhang et al. 1997, Rauch et 
al. 1997). In this paper we examine whether simulations like those of HKWM can also account 
for the observed Hell absorption, and we consider their implications for the physical state of the 
absorbing gas. 

The basic physics of intergalactic Hell absorption is elegantly described by Miralda-Escude 
(1993). The ambient background of UV radiation produced by quasars and (perhaps) young 
galaxies keeps the diffuse hydrogen of the Lya forest in a highly photoionized state, and it also 
ensures that most diffuse helium is singly or doubly ionized. However, the UV background 
intensity is lower at the ionization energy of Hell (4 Rydbergs, 228A) than it is at the HI ionization 
energy (1 Rydberg, 912 A), and Hell absorption can therefore be significant even in regions where 
the HI density is very low. Hell absorption may be the only practical tool for directly observing 
regions that lie significantly below the cosmic mean density, revealing gas whose hydrogen Lya 
absorption might be buried in noise or removed in the process of continuum fitting. According to 
gravitational instability models of structure formation, the voids between galaxies should harbor a 
substantial portion of the baryonic matter in the universe. A homogeneous intergalactic medium 
(IGM) would produce a uniform absorption trough in quasar spectra (Gunn & Peterson 1965). In 
a realistic gravitational instability model, the matter in underdense regions should instead produce 
a fluctuating continuum of absorption (Reisenegger & Miralda-Escude 1995). 

Once the hydrogen in the universe has been reionized, the neutral helium fraction is expected 
to be small except in high density, collapsed regions. At fixed density and temperature, the 
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fraction of Hell in highly photoionized gas is inversely proportional to the photoionization rate, 



where J(y) is the specific intensity of the background at frequency is, a^eii is the ionization cross 
section for Hell, and huueii = 4 Rydbergs. The HI fraction is inversely proportional to 



which is dominated by photons with energy between 1 and 4 Rydbergs. Measurements of the 
mean Hell and HI Lya opacities can constrain the spectral shape of the UV background through 
the ratio of these integrals, Tni/THeH, provided one has a model that specifies the density and 
temperature structure of the absorbing gas. The evolution of the UV background's intensity and 
shape can be tracked by the evolution of the mean opacities with redshift. The simulations provide 
a realistic model for the IGM with the physical detail needed to exploit this approach, albeit a 
model whose properties depend (as they should) on the parameters of the underlying cosmological 
scenario. 

For a given cosmological model, or indeed for any model IGM in which most of the 
absorption arises in highly photoionized gas, the predicted mean opacity depends on the parameter 
combination fif/r, where is the baryon density parameter (e.g., Miralda-Escude and Ostriker 
1992). This scaling assumes that gas temperatures (and hence recombination rates) are unaffected 
by changes in 0,b and T, an assumption that we will revisit in §3.2. The simulations allow us to 
ask whether a cosmological model is consistent with the observed mean opacities given constraints 
on Ofc from big bang nucleosynthesis (e.g., Walker et al. 1991) and on T from the observed quasar 
population (e.g., Haardt & Madau 1996) or the proximity effect (e.g., Bajtlik, Duncan, & Ostriker 
1988). The simulations also predict distribution functions for HI and Hell opacities and for the 
ratio of these opacities, which can be compared to observations in order to test the simulated IGM. 

The possiblity of detecting Hell absorption from diffuse intergalactic gas was one of the main 
scientific motivations for the Hopkins Ultraviolet Telescope (HUT), which first flew on the Astro- 1 
mission in 1990 (Davidsen, private communication). While the quasar absorption experiment could 
not be carried out during Astro-1, because of pointing problems, it was successfully performed 
during the Astro-2 mission in 1995, as described by Davidsen, Kriss & Zheng (1996, hereafter 
DKZ). The first detection of Hell absorption was in fact achieved before Astro-2 by Jakobsen et 
al. (1994, hereafter JBDGJP), who used the HST Faint Object Camera to observe the quasar 
Q0302-003 (z=3.28). They measured a clear drop in the received quasar continuum across the 
304A (rest-frame) edge, and they inferred a high Hell optical depth, THeii > 1.7 at 90% confidence. 
Because of the relatively low spectral resolution, the authors were unable to establish whether the 
absorption was caused mainly by material associated with individually identified HI lines or by 
a more diffuse component. JBDGJP, Giroux, Fardal, & Shull (1995), Madau & Meiksin (1995), 
and Songaila, Hu, & Cowie (1995) have explored the implications of these data for the shape of 
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the UV background spectrum assuming different analytic models of the absorbing medium. If 
the absorption is dominated by discrete lines, then the spectrum must be quite soft, but a harder 
spectrum is allowed if most Hell absorption arises in a diffuse, "Gunn-Peterson" background. We 
qualify this latter term with scare-quotes because the analytic models typically assume a uniform 
Gunn-Peterson effect, while the cosmological simulations predict a smoothly fluctuating, diffuse 
IGM that blurs the traditional distinction between Gunn-Peterson absorption and the Lya forest 
(HKWM; see also Miralda-Escude & Rees 1993; Cen et al. 1994; Miralda-Escude et al. 1996). 

Q0302-003 was recently reobserved by Hogan et al. (1996, hereafter HAR), using the Goddard 
High Resolution Spectrograph on HST. HAR confirm JBDGJP's detection of Hell absorption, 
but the higher spectral resolution of the GHRS observations reveals interesting new details. The 
transmission below the Hell break remains fairly high (~ 0.3) within about 4000 km s _1 of the 
quasar. HAR attribute the relatively low Hell fraction in this region to photoionization caused 
by Q0302-003 itself, with the ionization zone terminated by a high column density HI absorber. 
They find a low but significant level of residual flux at shorter wavelengths corresponding 
to THeii = 2.0^Q5 at 95% confidence. The upper limit on th c ii implies that helium remains 
predominantly doubly ionized even outside of the ionization zone produced by Q0302-003, 
presumably because of the ambient UV background. The upper limit depends on accurate 
background subtraction at the blue end of the spectrum, a challenging problem that is discussed 
in detail by HAR. 

HST cannot probe Hell absorption below z k, 3. DKZ took advantage of the shorter 
wavelength sensitivity of the HUT to measure Hell absorption in the spectrum of HS1700+64 
(z=2.743). They find a mean opacity in the redshift range 2.2 < z < 2.6 corresponding to 
THeii = 1-0 ± 0.07, with some evidence that the opacity increases as a function of redshift over this 
range. The observed wavelength interval is ~ 150A and the spectral resolution ~ 3A, so the HUT 
spectrum contains enough information to reveal structure in the residual flux. (The spectrum 
shown in DKZ is averaged in 10A bins.) The analysis in this paper will be aimed primarily at the 
HUT observations. The physical issues are rather different for Q0302-003 because of the important 
role of the observed quasar itself in ionizing the absorbing gas. We will therefore save a detailed 
comparison to the JBDGJP and HAR observations for a future paper. 

The numerical approach in this paper will be similar to that of HKWM, who use smoothed- 
particle hydrodynamics (SPH) simulations to study the HI Lya forest in a critical density, CDM 
universe. Studies of Hell (and HI) absorption in Eulerian-grid hydrodynamic simulations have 
been carried out by Zhang et al. (1995, 1996) and Miralda-Escude et al. (1996). Here we will study 
the physical state of the gas that produces Hell absorption in CDM-dominated, gravitational 
instability models of structure formation, relating it to and differentiating it from the gas that 
dominates HI absorption. We will use information on Hell and HI in the context of these models 
to study the required UV background spectrum, its evolution with redshift, and the implied 
baryon density. We will examine several variants of the CDM scenario, enabling us to see which 
features are generic within this cosmological picture and which can be used as diagnostics for 
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constraining cosmological parameters. We will also compare the simulation results to those of a 
simple model where all of the flux is absorbed in discrete lines. This sort of model is often used as 
a phenomenological description of the HI Lya forest, though it does not correspond physically to 
what happens in the simulations, where the absorbing structures are relatively diffuse and merge 
continuously into a fluctuating background. We obtain predictions that can be compared to future 
observational analyses that probe HI and Hell absorption along a common line of sight. 



2. Simulations 

We have used the N-body plus SPH code TreeSPH (Hernquist & Katz 1989; Katz, Weinberg, 
& Hernquist 1996, hereafter KWH) to simulate three different CDM-dominated cosmological 
models, the parameters of which are listed in columns 2 — 5 of Table 1. The first is a "standard" 
CDM (SCDM) universe, with n = 1, h = 0.5 (where h = H /W0 km s" 1 Mpc" 1 ), and n b = 0.05. 
The power spectrum is normalized so that the rms amplitude of mass fluctuations in 8 h~ 1 Mpc 
spheres, linearly extrapolated to z = 0, is eg = 0.7. This normalization is consistent with that 
advocated by White, Efstathiou & Frenk (1993) to match the observed masses of rich galaxy 
clusters, but it is inconsistent with the normalization implied by the COBE-DMR experiment. 
Our second model is identical to the first except that as = 1.2. This higher amplitude is consistent 
with the 4- year COBE data (Bennett et al. 1996), and we therefore label the model CCDM. 
The third model, OCDM, assumes an open universe with fio = 0.4, h = 0.65 and f2& = 0.03. 
The transfer functions used are those of Efstathiou, Bond & White 1992. The shape parameter 
T = 0.234 was used for the OCDM model, which is also COBE-normalized (erg = 0.75, Gorski et 
al. 1995). The baryon fraction for these models, = 0.0125/i~ 2 , was chosen based on the big bang 
nucleosynthesis analysis of Walker et al. (1991), who deduce primordial abundances of D, 3 He, 
4 He, and 7 Li from local observations using models for chemical processing of elements in stars. 
Measurement of the deuterium abundance in high redshift Lyman limit systems offers a more 
direct route to determining O^/i 2 . The first applications of this method have so far given results 
which favor Q^h 2 being a factor of ~ 2 smaller (Carswell et al. 1994; Songaila et al. 1994; Rugers 
& Hogan 1996ab, but see also Hogan 1997) or a factor of ~ 2 larger (Tytler, Fan Sz Buries 1996, 
Tytler, Buries & Kirkman 1996, but see also Songaila, Wampler & Cowie 1996). The observational 
situation is thus uncertain at present, though we shall see that the high tl^h 2 results are favored 
in the CDM models. 

A periodic cubic volume of comoving side length 11.111 /i _1 Mpc was simulated for each 
model, each simulation having identical random phases and being evolved to z = 2. We analyzed 
outputs at z = 3.5, 3.0, 2.67, 2.33 and 2, encompassing the range of existing Hell absorption 
observations. Each simulation was run with 64 3 collisionless dark matter particles and 64 3 gas 
particles. A spatially uniform photoionizing radiation field was imposed, and radiative cooling 
and heating rates were calculated assuming photoionization equilibrium and optically thin gas, as 
discussed in KWH. The spectral shape of the background radiation field and its evolution with 
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redshift were taken from the work of Haardt & Madau (1996, hereafter HM). This background 
was calculated by HM from a self-consistent treatment of the absorption and re-emission of light 
from observed quasars by the Lya forest. 

The implementation of photoionization in TreeSPH requires values of the photoionization 
rates and photoionization heating rates for HI, Hel, and Hell as a function of redshift (see Katz, 
Weinberg, & Hernquist 1996). We compute these parameters from the HM spectrum at intervals 
Az = 0.05, and for the simulations we use fitting formulae that match these values to within 
~ 10%. The parameters that are particularly relevant for this paper are the photoionization rates 
rni and THeii, for which we adopt 

HM r ffl = 1.115 x 10" 12 exp [-0.57565(z - 2.5) 2 j s _1 , (3) 

HMr H eii = 1-088 x 10~ 14 exp [-0.57565(z - 2.5) 2 j P(z) s' 1 , (4) 

where the factor 

P(z) = 1 + 0.125(2 - 2.5) + 0.0825(z - 2.5) 2 (5) 

accounts for the slight difference in the relative evolution of Thi and Tneii- At z ~ 2 — 4, the values 
implied by equations (|3|) and (||) are about 30% lower than those shown in figure 6 of HM because 
we fit HM's qo = 0.5 results rather than their q$ = 0.1 results (see §5.13 of HM). We set the UV 
background to zero at z > 6, since observations suggest that the population of quasars is already 
declining rapidly between z = 4 and z = 5. As shown by Hui & Gnedin (1997), the thermal state 
of the IGM at z ^ 3 is insensitive to the epoch of reionization provided that it occurs at z ^ 5. If 
HI or Hell reionization occurs at z < 5, the IGM could be somewhat hotter than our simulations 
predict, a point that we will return to at the end of §3.4. 

The HKWM results suggested that a simulation with the HM background and = 0.05 
would underproduce the observed HI Lya opacity, and we therefore reduced the amplitude of the 
HM background by a factor of two (i.e. to half the values implied by equations 3 and 4) before 
evolving the simulations. We can recover the impact of different UV background intensities at the 
analysis stage, regardless of the specific value adopted during dynamical evolution, as discussed 
below in §3.1. The SCDM simulation of this paper is identical to the simulation analyzed by 



Table 1: Parameters of cosmological models. The conventions are defined in §2, except for C 
which is defined in equation (7). 

model r?o h fi& o"s ChoII (z = 2 — 3.3) Cui(z = 2) Cui(z = 3) 



SCDM 1 0.5 0.05 0.7 
CCDM 1 0.5 0.05 1.2 
OCDM 0.4 0.65 0.03 0.75 



3.0 
8.0 
3.0 



2.6 
5.7 
2.4 



2.8 
6.7 
2.8 
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HKWM except for the UV background spectrum (HKWM used a v~~ l power-law) and the inclusion 
of star formation (using the prescription of KWH). Star formation and the associated feedback 
only influence the simulation results in high density regions that are not the focus of this paper. 

The UV background determines the relative fractions of different ionic species at specified 
density and temperature. Figure 1 shows the fractions of HI, Hel, and Hell as a function of 
temperature for (total) hydrogen densities uh = 3.95 x 10~ 6 cm' 3 (left panel) and 3.95 x 10~ 5 cm' 3 
(right panel). These are, respectively, the mean density and an overdensity of 10 at z = 2.33 for 
n b h 2 = 0.0125. We will focus much of our analysis on this redshift because it is the simulation 
output closest to the middle of DKZ's observed redshift range. Abundances are computed by 
requiring balance between the production and destruction rates for each species, as described 
in KWH. Thin lines show abundances for no UV background — these collisional equilibrium 
abundances depend on temperature alone and are thus the same in the two panels. Thick lines 
show abundances for the HM background, with intensity reduced by a factor of three (relative to 
HM), the value required for our SCDM and OCDM models to reproduce the observed Hell opacity 
(see §3.1). Photoionization dominates completely over collisional ionization at low temperatures, 
where the HI and Hell fractions fall as T' 0,7 because of the temperature dependence of the 
recombination rates. Collisional ionization becomes significant at high temperatures. Raising the 
density increases the recombination rates and, therefore, the HI, Hel, and Hell fractions, though 
these never exceed their collisional equilibrium values. Lowering the UV background by a constant 
factor would have the same effect on these plots as raising the density by the same factor. 

The shape of the background spectrum determines the relative fractions of HI and Hell 
through the ratio of photoionization rates, Thi/ThcH- HM assume a u^ 1 - 5 power-law for the 
intrinsic UV spectrum of their quasar sources, but the ambient background that they compute 
is substantially softer because of reprocessing by the Lya forest. In particular, the spectrum at 
z ~ 2.5 drops by a factor ~ 15 at hv ~ 4 Rydbergs because of Hell absorption. If we adopted 
a pure z/ -1,5 power-law with intensity chosen to produce the HI fractions in Figure 1, the Hell 
fractions at low temperatures would be reduced by a factor of 3.3. With the v~ l power-law of 
HKWM, the reduction would be a factor of 6.5. It is worth noting, however, that these changes 
in equilibrium abundances have no dynamical impact on the simulations and minimal impact on 
the ability of the gas to cool in collapsed objects. The fraction of gas that cools and condenses 
into galaxies is nearly identical in simulations with the HM background and ai/ _1 background 
(Weinberg, Hernquist, & Katz 1996). 



3. Results 

The Lya optical depth produced by a uniform medium with neutral hydrogen or singly 
ionized helium density n is (Gunn & Peterson 1965; Miralda-Escude 1993) 

r = —f\H-\z)n 1 (6) 
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Fig. 1. — Fractions of HI (solid lines), Hel (dotted lines), and Hell (dashed lines) as a function 
of gas temperature. The assumed gas densities are njj = 3.95 x 10~ 6 cm~ 3 in the left panel (the 
mean density at z = 2.33 for Vt^h? = 0.0125) and a factor of ten higher in the right panel. Thin 
lines show collisional equilibrium abundances. Thick lines show abundances in photoionization 
equilibrium with a UV background that has the spectral shape computed by HM but an intensity 
lower by factor of three. 
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where H(z) is the Hubble constant at redshift z, f = 0.416 is the Lya oscillator strength, and A 
is the transition wavelength (1216A and 304L4 for HI and Hell, respectively). In the simulations, 
the intergalactic gas is not uniform, so the optical depth varies as a function of wavelength and 
angular position. Figure 2 shows HI and Hell "absorption maps" of a slice through the center 
of the SCDM simulation at z = 2.33. (the velocity width of the full cube is 2029 km s _1 at 
this redshift). The faintest grey levels correspond to r ~ 0.05, a value commonly adopted as a 
threshold for line identification in high signal-to-noise spectra. Absorption saturates for r ^ 3 
(e-^0.05). 

As expected, these maps show stronger Hell absorption than HI absorption. Hell produces a 
substantial optical depth in regions that are nearly devoid of HI absorption, so Hell measurements 
can probe low density structure that is virtually undetectable by other means. Conversely, most 
regions with a significant HI optical depth produce saturated Hell absorption. The sharp edges 
visible in some of the saturated Hell features and the "blobbiness" of the weak Hell absorption 
are probably artifacts of representing the gas distribution with a discrete set of particles. 

For quantitative analyses, we create simulated absorption spectra along random lines of sight 
through the simulation cube. Knowing the temperature and density of each SPH particle, we 
can compute the associated fractions of HI, Hel, and Hell given our assumed UV background. 
From these we compute the optical depths as a function of frequency by performing a line integral 
through the smoothing kernels of all SPH particles whose kernels intersect the line of sight, taking 
into account Hubble flow, peculiar velocities, and the thermal broadening appropriate to each of 
the absorbing species in order to transform from physical space to frequency space. Details of this 
procedure are described in HKWM. 

3.1. Mean absorption and evolution with redshift 

Figure 3a shows the redshift evolution of the effective mean optical depth for Hell absorption, 
fueii = — log e (F}, where F is the transmitted flux (with F = 1 corresponding to complete 
transmission). The average is performed over 100 lines of sight at each of the five output times. 
Because of the non-linear dependence of flux on optical depth, THeii is not the mean value of 
THeii along the spectrum, but while (F) can be measured from a spectrum of imperfect resolution 
and finite signal-to-noise ratio, the mean value of th c ii cannot. Thin lines show results for the 
three models with the UV background employed during dynamical evolution of the simulation, 
i.e., the HM background divided by a factor of two. The intensity of this UV background at 
energies responsible for ionizing Hell is approximately constant over the redshift range plotted 
here, though it does start to decline smoothly above z ~ 3.5. The mean optical depth decreases 
towards lower redshift primarily because the expansion of the universe reduces the Hell fraction, 
since it lowers the physical density of the absorbing medium and hence the recombination rate per 
Helll ion. Cosmic expansion also spreads the atoms in a given comoving region over a larger range 
in frequency. These are precisely the effects that determine the evolution of the Gunn-Peterson 
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Fig. 2. — Optical depth to HI (left) and Hell (right) absorption in the SCDM simulation at 
z = 2.33. The optical depth is averaged over a slice of velocity width 60 km s _1 through the 
middle of the simulation cube. The scale is logarithmic in r, as indicated, with the faintest grey 
levels corresponding to r ~ 0.05. The intensity of the UV background has been scaled so that 
the simulation produces the mean Hell opacity measured by DKZ (see discussion in §3.1). Hell 
absorption is prominent in low density regions that are nearly transparent in HI. 
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optical depth produced by a uniform IGM with a constant UV background. The first effect, which 
dominates here, would not apply to gas in collapsed, physically stable structures. While Figure 2 
shows that the absorption is by no means uniform, the rapid evolution of Tjj e n is a strong hint 
that the absorbing gas resides in low density structures that are still expanding with the cosmic 
background. 

The solid triangle in Figure 3a shows HAR's observational estimate, rneii = 2.0^^ (2a 
errors) at z = 3.3, measured outside the ionization zone produced by the observed quasar. The 
wavelength range contributing to this point is about 40A Solid circles with error bars show fHeii 
at lower redshifts from the DKZ spectrum. Each point corresponds to a single 10A bin from 
DKZ's figure 1 — we measured the extrapolated quasar continuum and observed flux from this 
figure and divided them to obtain F = exp(— rneii)- The error bars are la and are based on the 
noise vector in DKZ's figure. The fluctuations from point to point are larger and more coherent 
than expected from observational errors alone. They are presumably caused by large scale density 
variations, or possibly by inhomogeneity in the background radiation field. 

Comparison between the data points and the thin lines in Figure 3a shows that the CDM 
models produce roughly the correct trend of rneii with redshift but have an overall level of 
absorption that is too low. The intensity of the ionizing background used during the dynamical 
evolution of the simulations is already a factor of two lower than that advocated by HM, but a 
further reduction is necessary to match the observed mean absorption. To quantitatively compare 
the models, we calculate the factor 

n HmFhcII 

C/Hell = ^ , (7 J 

reql Hell 

where reqFHcii is the photoionization rate required to match the observed Hell optical depth 
and HMTHeii is the photoionization rate predicted by the HM background (equations 3 and 4). 
Changing the intensity of the UV background within a broad range does not significantly alter the 
temperatures or spatial distribution of the absorbing gas, so we are able to estimate reqFHeii by 
rescaling the background after the simulation has been run but before creating simulated spectra 
(see Miralda-Escude et al., in preparation, for numerical tests of this scaling). In fact, because 
photoionization completely dominates over collisional ionization in the absorbing regions, we could 
obtain virtually identical results by the still simpler procedure of rescaling optical depths in the 
extracted spectra. 

The Hell absorption value we match is the quoted value of rneii = 1-0 =t 0.07 at z = 2.4 from 
DKZ. We find Cneii by an iterative interpolation search, dividing the HM ionization parameters 
by different factors until we obtain the correct mean absorption. As we do not have a simulation 
output at z = 2.4, we linearly interpolate rneii between z = 2.33 and z = 2.67. The values of 
Cfleii for the three cosmological models are listed in column 6 of Table 1. These factors are equal 
for the two low amplitude models, OCDM and SCDM. The required value for the CCDM model 
is much higher. This model tends to produce less absorption at fixed V because (as we will show 
in detail later) most of the flux decrement arises in large volumes that are near or below the mean 



-12- 




Fig. 3. — (a) Redshift evolution of THeii; the effective mean optical depth of Hell absorption. The 
horizontal axis (redshift) is scaled logarithmically in (1 + z), so that power-law relations between r 
and (1 + z) would appear as straight lines. Thin lines with open squares show the results for the 
three CDM simulations calculated with the ionizing background intensity used during dynamical 
evolution. Thick lines with filled squares show THeii after the intensity of the background is reduced 
by a factor of 1.5 for SCDM and OCDM and by 4.0 for CCDM. With this lower intensity, the 
simulations reproduce the mean optical depth measured by DKZ, THeii = 1-0 at z = 2.4. Points 
from the individual lOA bins of the DKZ spectrum are shown as solid dots with la error bars. 
The triangle shows the measurement and 2a error bar of HAR, th c ii = 2.01^ at z « 3.3. (b) 
Same as (a), but for HI absorption. Thin lines again correspond to the UV background used during 
dynamical evolution and thick lines to the background reduced by a factor of 1.5 (SCDM, OCDM) 
or 4.0 (CCDM). The dotted line is PRS's fit to the observed evolution of fm, which was calculated 
from data in the redshift range indicated by the hatched region. The vertical extent of this region 
illustrates the approximate la confidence interval for the fit. The solid points are the measurements 
of mean optical depth from ZL for their "sample 2" quasars. 
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density. At the higher mass fluctuation amplitude of CCDM, more of the gas has flowed out of 
these void-like regions into higher density zones and collapsed objects; the absorption from these 
regions is already saturated, so adding more gas to them does not increase the mean decrement. 

The heavy lines in Figure 3a show THeii with the rescaled ionization parameters. A single 
value of Cneii allows the models to fit the observed THeii results at all redshifts from z = 2 to 
z = 3.3, to within the scatter of the data points, indicating that the evolution predicted by the 
simulations is consistent with the observations if the UV background evolves as predicted by HM. 
Although Tneii falls by only 20% between z = 2.4 and z = 3.3, the models reproduce the factor of 
two difference between the DKZ and HAR measurements of THeii because of the cosmic expansion 
effects discussed earlier. 

Figure 3b shows the effective mean optical depth of HI absorption, thi = — log e (FHi), for 
the three CDM models. Thin lines correspond to the UV background intensity adopted in the 
simulations, i.e., half of the HM background intensity at the corresponding redshift. Thick lines 
show thi with the UV background intensity rescaled in order to match the observed rneiij i-e., the 
HM background (equations 3 and 4) divided at all redshifts by the factor ChcII listed in Table 1. 

The straight dotted line in Figure 3b shows the power-law fit thi = 0.0037(1 + z) 3 ' 46 found 
by Press, Rybicki, & Schneider (1993, hereafter PRS) in their analysis of the mean flux decrement 
in a sample of 29 high-redshift quasars. The hatched box shows the redshift extent of points 
used to compute this fit and the la statistical uncertainty in the fit. PRS computed the mean 
decrement by extrapolating the quasar continuum from the region redward of Lya emission into 
the Lya forest region. Zuo &: Lu (1993, hereafter ZL) estimated the mean decrement from higher 
resolution spectra by directly fitting the continuum to the regions of lowest absorption in the 
forest. We indicate their results by the filled circles in Figure 3b. Clearly the two observational 
determinations disagree by far more than their statistical uncertainties. A recent analysis of Keck 
HIRES spectra, using an approach similar to ZL's but data of higher resolution and signal-to-noise 
ratio, yields mean decrements that are much closer to the PRS values than to ZL's (Rauch et al. 
1997). We will therefore proceed on the assumption that the PRS determination is accurate, but 
the systematic uncertainty in existing estimates of thi is worth keeping in mind, as a major change 
to these estimates would affect our conclusions about the shape of the UV background spectrum. 

Once the background intensity has been divided by the factor Ch c ii, the effective HI optical 
depths agree fairly well with the PRS determination, though they tend to rise above it at the low 
and high ends of our redshift range. To better quantify this agreement, we list in columns 7 and 
8 of Table 1 the scaling factors Chi (defined analogously to Ch c ii) by which the HM background 
must be divided in order that the simulation match the PRS optical depths at z = 2 and z = 3. 
These factors match the corresponding values of Cneii to 20% or better in most cases (30% for 
CCDM at z = 2), indicating that the cosmological simulations and the HM spectral shape are, 
taken together, consistent with the joint observational constraints of DKZ and PRS. The HI 
scaling factors at z = 2 and z = 3 are also similar, indicating that the simulations reproduce the 
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PRS evolution law over this redshift range if the UV background evolves as predicted by HM. 
If this analysis is extended to z ~ 4, however, the simulations require a roughly constant Thi 
(HKWM; Rauch et al. 1997), while the HM model predicts a substantial drop in Thi towards high 
redshift because of the declining number density of quasar sources. 

The uncertainty on the correct value of Cneii is mainly statistical, and is dominated by the 
fact that an observational measurement of Tjjeii is only available from one QSO. The theoretical 
estimates at z=2.33 come from 200 lines of sight through the simulation box, which together cover 
a redshift patch equal to ~ 13 times the useful length of the DKZ spectrum. We can estimate 
the error in Cneii by picking groups of simulated spectra with the same total length as DKZ from 
our ensemble of 200. We then calculate the CHeii required to match the observational THeii for 
each set of spectra. From the spread of Cneii values we estimate the la uncertainties on Ch c ii to 
be +10%, —25%. The error in Chi is dominated by systematic uncertainties in the observational 
determination of thi, probably associated with continuum fitting (see Rauch et al. 1997). 

JBDGJP suggested that the high Hell optical depth towards Q0302-003 might arise because 
the Hell "Stromgren spheres" around quasars had not overlapped by z = 3.3, leaving most of 
the universe optically thick to Hell ionizing photons. Supporting evidence for this scenario comes 
from a rapid change at z ~ 3.1 in the SilV/CIV ratios measured in Lya absorbers, which suggests 
a change in the shape of the UV background spectrum that could correspond to percolation of 
quasars' Hell ionization zones (Songaila & Cowie 1996, but see also Hellsten et al 1997). Our 
results show that, in the CDM models, no major change in the spectral shape is needed to explain 
the existing Hell data. Indeed, if HAR's upper limit on THeii is taken at face value, then the 
background spectrum cannot be much softer than the HM spectrum at z = 3.3, at least along this 
line of sight. If HAR's residual flux is an artifact of imperfect background subtraction, then our 
models could also accommodate Hell reionization at z < 3.3. 

3.2. Implications for 

All of our models require Chi ^ 2.5 to match the PRS determination of thi and Cneii ^ 3 to 
match the DKZ determination of THeii • If we had obtained values of Chi and Cneii smaller than 
unity from the simulations, we could accommodate them easily by appealing to additional UV 
sources not considered by HM, e.g., star-forming galaxies or faint AGNs. However, since the HM 
background is based on the observed population of quasar sources (with a modest extrapolation for 
quasars below existing survey detection limits), it is difficult to see how the true background could 
be lower than the HM background by such large factors. Reductions of this magnitude would also 
make the background intensity inconsistent with estimates from the proximity effect (Giallongo et 
al. 1996 and references therein), though these are subject to significant systematic uncertainties. 
Even allowing for the imperfect resolution of the simulations and plausible uncertainties in the 
background intensity, it seems that these cosmological models for the Lya forest are at best 
marginally compatible with the PRS opacity measurements if f^/i 2 = 0.0125. 
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The alternative to lowering the UV background is to raise the mean baryon density. The 
Lya optical depth is proportional to the number density of absorbing atoms (see equation ||), 
which for highly photoionized gas is proportional to the square of the gas density divided by the 
photoionization rate. The second power of density arises because the the recombination rates per 
HII or Helll ion are themselves proportional to the density. If the distribution of oue?densities 
(Pb/Pb) and gas temperatures in the IGM is unchanged by altering Q&, then the optical depth 
r at a specified redshift along a line of sight is proportional to £1%/T. Note that the effective 
mean optical depth r is not simply proportional to f^/r because of the non-linear nature of flux 
averaging, but it is still the case that raising £lb by a factor C 1//2 should have the same effect as 
lowering T by a factor C. The assumption that the absorbing gas is highly photoionized breaks 
down in regions that are collisionally ionized (e.g., hot gas in virialized halos) or predominantly 
neutral (e.g, damped Lya systems), but these are too rare to make much contribution to r. 

We have completed one simulation of the SCDM model with Qf, = 0.125 instead of 0.05, and 
we find that the scaling of r with £lb is weaker than the above argument would suggest, roughly 
r oc f^' 7 instead of r oc The reason is that raising Qj, also raises the gas temperature in the 
low and moderate density regions (pb/pb ^ 10) that produce most of the absorption, because 
increasing the HI and Hell fractions allows a given volume of gas to absorb energy from the 
photoionizing background at a higher rate. Since the recombination rates decline as T~ - 7 in the 
relevant temperature regime, the HI and Hell fractions do not rise by the full factor when f2& 
is increased. The physics of the Qb scaling does not depend on the cosmological scenario, and we 
therefore expect the result derived from our pair of SCDM simulations to hold more generally 
(but see the discussion of reionization effects at the end of §3.4). One might think that reducing T 
at fixed Ob would also alter the gas temperatures, but it does not, because the increase in HI and 
Hell fractions is exactly compensated by the smaller rate of photoionizations per ion. The scaling 
r oc T" 1 is therefore preserved unless the shape of the ionizing background, which determines the 
mean residual energy per photoelectron, is altered. 

With these scalings in mind, we can relate the values of Chi and Cneii listed in Table 1 to 
the combination of T, Qf,, and h that is required for the simulation to match the observed mean 
opacity: 

/ \ 1 7 

/£hm\ / &bh \ fh sim 



c -{t^-j{omk) brj' (8) 

where /i s j m is the value of h adopted in the simulation. The h dependences arise because the 
mean gas density is proportional to h? at fixed Qb and because the optical depth at fixed HI or 
Hell density is inversely proportional to the Hubble constant (equation ||). If we assume that 
r true = Thm and h = h s i m , then matching the PRS values of thi requires Q,bh 2 ~ 0.023 for SCDM 
and OCDM and Qbh 2 ~ 0.038 for CCDM. Matching the DKZ measurement of rneii requires a 
similar baryon density, though the Hell opacity on its own gives a less compelling argument for 
high Ofe because the lower limit on Tneii is less secure than the lower limit on Thi- Assuming 
standard big bang nucleosynthesis, a density £lbh 2 = 0.023 is in excellent agreement with Tytler 
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et al.'s (1996) estimate of the primordial deuterium abundance, but it is inconsistent with the 
much higher deuterium abundances estimated by Carswell et al. (1994), Songaila et al. (1994), and 
Rugers & Hogan (1996ab). We examine the baryon density required by cosmological simulations 
of the Lya forest more thoroughly in Rauch et al. (1997), which includes a new determination 
of Tm{z) from high resolution spectra and a discussion of the lower limit on Thi- The analytic 
arguments presented in Weinberg et al. (1997) show that the lower limits on derived from this 
method depend only on very general properties of the "cosmological" picture of the Lya forest, 
and are unlikely to be weakened substantially by changes in the adopted cosmological model or 
the numerical resolution of the simulations. 

3.3. Simulated Hell spectra 

The top panels of Figure 4 show examples of HI (solid lines), Hel (dotted lines), and Hell 
(dashed lines) absorption along two randomly selected lines of sight through the SCDM simulation, 
at z = 2.33. We have scaled the HM ionizing background by the factor Cneii listed in Table 1, 
so that the mean optical depth matches the DKZ observation. The transmission e _T is plotted 
against line-of-sight velocity. The corresponding baryon density (in units of the mean baryon 
density) is plotted below each spectrum in the second panel, with the solid line showing the 
redshift space density and the dotted line the real space density (i.e., the density computed with 
peculiar velocities and thermal broadening set to zero). The HI optical depth is well correlated 
with the redshift space density. Features in the redshift space density field are usually offset from 
those in the real space field because of peculiar motions, and the high density peak in the left hand 
spectrum is greatly broadened in redshift space because of infall. Hel absorption is non-negligible 
only in the highest density region of the second spectrum, where the high recombination rate 
increases the relative fraction of neutral helium. Hell absorption, on the other hand, is quite 
strong, and most of it arises in regions where the HI optical depth is quite low. The ratio t Hell / Thi 
(shown by the solid line in the third panel) is about a factor of eight over most of the spectrum. 
The HI and Hell flux decrements, 1 — e~ T , are nearly equal when the HI optical depth is high, but 
when thi is small the Hell flux decrement is eight times higher. Hell absorption therefore probes 
regions of lower density than HI absorption, as already seen in Figure 2. 

The optical depth ratio THeii/THi is high and approximately constant, as anticipated by 
Miralda-Escude (1993). Variations arise when collisional ionization or thermal broadening of 
the spectrum become important — these processes affect Hell and HI differently because of the 
differences in ionization potential and atomic mass, respectively. The dotted line in the third 
panel of Figure 4 shows the optical depth ratio calculated from spectra along the same lines of 
sight with no thermal broadening applied. The variations seen in the solid line largely disappear, 
indicating that thermal broadening is their primary cause, at least along these two lines of sight. 

For a more quantitative view of the optical depth ratio, we examine the joint distribution of 
THcii and thi in 200 spectra extracted along random lines of sight through the SCDM simulation 
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Fig. 4. — Two examples of spectra taken along random lines of sight through the simulation box of 

the SCDM model, at a redshift z = 2.33. The top panel shows transmitted flux against line-of-sight 

velocity for HI (solid line), Hel (dotted line), and Hell (dashed line). The second panel shows the 

density of baryonic matter, in units of the mean baryon density, in redshift space (solid line) and 

real space (dotted line). The third panel shows the ratio THeii/Tn computed with (solid line) and 

without (dotted line) thermal broadening of the absorption. The bottom panel shows the ratio 

A = men/pi' 6 computed in redshift space (solid line) and real space (dotted line). 
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Fig. 5. — The joint distribution of THeii and tri in the SCDM model at z = 2.33, computed from 
200 spectra, each with 1000 pixels. The logarithmic grey scale shows the fraction of pixels in bins 
of Athi = 0.12, AthcII = 1.2. Lines show the percentile distribution of TReii in each bin of tri- 
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Fig. 6. — The joint distribution of optical depth and pi> (in units of the mean baryon density) for 
the SCDM model at z = 2.33, in a format similar to Fig. 5. (a) Hell in redshift space, (b) HI in 
redshift space, (c) Hell in real space, (d) HI in real space. 
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at z = 2.33. The logarithmic grey scale in Figure 5 indicates the distribution of pixels in the 
Thi — THeii plane. There are 2 x 10 5 pixels in total, 1000 in each of the 200 spectra. The labeled 
lines in Figure 5 show the percentile ranges of THeii in bins of Athi = 0.12, i.e., of the pixels 
that have a given value of thi, 1% have THeii below the 1% line, 5% below the 5% line, and so 
forth. Most pixels lie along a well defined ridge at THdi = 8thi, tracked by the median 7Heii line. 
Collisional ionization in shock heated regions raises THeii / Thi because it suppresses HI absorption 
more strongly than Hell absorption. Thermal broadening on the edges of high density regions 
tends to reduce THcIi/thi by spreading HI absorption into the lower density surroundings. Both 
effects can be seen in Figure 5, but the scatter below the median relation is broader than the 
scatter above it, confirming the anecdotal evidence of Figure 4 that thermal broadening is the 
dominant cause of variations in THeii/ Thi- 

Most of the absorbing gas in Figure 4 has relative density p/p between 0.1 and 10. In 
these simulations, the gas in this density regime typically follows a tight and simple relation 
between temperature and density, approximately T cx p 0M (see figure 3b of Weinberg, Hernquist 
& Katz 1997). This relation arises because the gas is cooled by adiabatic cooling and heated by 
photoionization, at a rate that depends on the density. Changing the gas temperature alters the 
heating rate, and gas that lies off this relation evolves towards it on a Hubble timescale (for a 
more detailed discussion see Hui & Gnedin 1996). The relation breaks down in collapsed regions, 
where shock heating and radiative cooling become important. The optical depth to HI or Hell 
Lya absorption is proportional to the number density of HI or Hell atoms at the corresponding 
line-of-sight velocity. These are proportional to the gas density multiplied by the recombination 
rate, which is in turn proportional to pT~ ' 7 for temperatures up to several x 10 4 K. Gas that lies 
on the T cx p ' 6 temperature-density relation should therefore satisfy th c ii /9 2 (p ' 6 ) -0 ' 7 cx p 1 ' 6 . 
The solid line in the bottom panel of Figure 4 plots the ratio A = THeii/p 1 ' 6 against velocity, 
and it is indeed approximately constant over most of the spectrum. The variations are caused 
primarily by peculiar velocity distortions of the redshift space density, since it is the real space 
density that is directly correlated with the temperature. In real space (dotted line) the variations 
are much smaller, with the one major departure occurring in the high density region of the left 
hand spectrum, where shock heating drives the gas off the simple temperature-density relation. 
Figure 6a shows the joint distribution of THeii and p, in redshift space, for the 200 spectra 
examined in Figure 5. Most pixels lie on the ridge THeii ~ 3.5(p/p) L6 , though there is a scatter 
towards lower rneii at higher densities because of shock heating. The joint distribution of thi and 
p (Figure 6b) is similar, except for a factor of eight offset. The factor of two interquartile scatter 
at low densities in these plots is caused predominantly by peculiar velocity effects, as one can see 
from the corresponding real space joint distributions (Figures 6c and 6d). 

To a first approximation, one can thus regard a Hell or HI absorption spectrum as a map of 
the gas density field along the line of sight, albeit a map that is non-linear (an exponential of a 
power-law) and distorted by peculiar velocities. Figures 4-6 use the SCDM model for illustration, 
but the qualitative physical picture is similar in all three of the cosmological scenarios that we 
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consider. This picture can be contrasted with a traditional phenomenological description of the 
Lya forest as a collection of discrete absorbing clouds, each producing a Voigt-profile line fully 
characterized by a redshift, an HI column density, and a 6-parameter (velocity width). One can 
compute the Hell absorption produced by "line blanketing" in such a model by assuming a ratio 
nueii/nm and either thermal broadening (in which case the Hell 6-parameters are a factor of two 
smaller than the HI 6-parameters) or "turbulent" broadening (in which case the 6-parameters are 
equal). For a heavily saturated but undamped line, the equivalent width is proportional to the 
6-parameter. 

In order to compare the simulation results to this sort of discrete cloud model, we have used 
a program written by J. Miralda-Escude to generate artificial spectra that are superpositions 
of randomly distributed, Voigt-profile lines. We draw HI column densities from a power-law 
distribution, dN/dNm oc N^' 5 , with a lower cutoff at iVHi,min = 10 12 cm~ 2 , and 6-parameters 
from a Gaussian distribution with a mean of 28 km s _1 and a dispersion of 10 km s _1 , truncated 
below 6 m i n = 18 km s _1 . These parameters, based on Hu et al. (1995), are similar to those used 
by Songaila et al. (1995) in their modeling of the JBDGJP observation, though we have lowered 
A^Hi,min from 2 x 10 12 cm~ 2 to 10 12 cm~ 2 . We choose the mean number of lines per unit redshift 
in order to match the PRS determination of thi at z = 2.33. We then generate two sets of 
corresponding Hell spectra, one for pure thermal broadening, one for pure turbulent broadening, 
choosing the nneii/^Hi ratio in each case to reproduce the DKZ value of th c ii at z = 2.33. 

Figure 7 shows two examples of the line model spectra. Comparing to Figure 4, we see that 
while the HI spectra look qualitatively similar to those produced by the cosmological simulations, 
the Hell spectra look quite different — they are sharply corrugated, with many saturated regions 
in each spectrum. The simulations reproduce the observed values of Thi and thcii at z = 2.33 
simultaneously if Thi/FroU ~ 100, as implied by the HM background spectrum. Our line models 
reproduce thi and THeii by construction, but the required UV background spectrum is much 
softer, with Thi/PhoU ~ 200 for the turbulently broadened model and ~ 2500 for the thermally 
broadened model. If the minimum column density is pushed far below 10 12 cm~ 2 , then the 
qualitative properties of the line model become closer to those of the simulations, since the weak 
lines overlap to produce a fluctuating background that gives rise to much of the Hell absorption. 
Whether the underlying physical picture approaches that of the simulations depends on how one 
envisions the absorbers themselves. In a discrete cloud model, the absorption arises in physically 
distinct objects whose wings overlap in frequency space because of line broadening, but in the 
cosmological simulations the absorption arises in a smoothly fluctuating, continuous IGM. 

3.4. Statistical analysis of transmission and optical depth 

We now turn to statistical measures that quantify properties of the absorbing gas in the 
various models that we have introduced. These statistical predictions can be used to test and 
differentiate these models. 
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Fig. 7. — Two spectra generated from a model of randomly distributed, Voigt-profile lines. Dashed 
lines in the top panels show transmitted flux in the HI Lya spectrum. Solid and dotted lines 
show Hell transmission assuming thermal broadening and turbulent broadening, respectively. The 
spectra are normalized to produce the observed values of thi and THeii • The bottom panels show 
the ratio THeiiAlu- This is constant for the turbulent broadening model (dotted line) because each 
line's HI and Hell 6-parameters are identical. In the thermal broadening model (solid line), the 
ratio varies strongly because of the factor of 2 difference in 6-parameters. 
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Fig. 8. — Mean absorption as a function of gas overdensity, at z = 2.33. D = (1 — F) is the mean 
flux decrement in the spectrum after the contribution of gas with density below a threshold value of 
Pb (in units of the mean baryon density) is eliminated. Thick lines show results for Hell, thin lines 
for HI. The SCDM, CCDM, and OCDM models are represented by solid, dashed, and dot-dashed 
lines, respectively. 
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Figure 8 demonstrates a point suggested qualitatively earlier: most of the Hell absorption in 
the cosmological simulations arises in low density regions. We compute the mean flux decrement, 
D = (1 — F), after setting to zero the absorption caused by gas with baryon density below a 
specified threshold. The flux decrement is plotted as function of the pb threshold (in units of the 
mean baryon density). Regions with overdensity pb/ pb < 2 account for roughly 80% of the Hell 
flux decrement. This density regime should be amenable to approximate analytic treatments - 
though not necessarily to linear perturbation theory per se — so analytic methods like that of 
Reisenegger &; Miralda-Escude (1995; see numerical tests in Miralda-Escude et al. 1996) should 
provide useful guides to the predictions of Hell absorption in CDM-like models. We see from 
Figure 8 that the SCDM and OCDM models, which have similar mass fluctuation amplitudes, 
produce their Hell absorption in almost identical density regimes. In the higher amplitude, CCDM 
model, noticeably less absorption arises in moderately underdense regions. As discussed in §3.1, 
more material in the CCDM model has flowed out of "voids" into relatively dense objects, which 
produce saturated Hell absorption. If the intensity of the UV background were held constant, this 
model would produce substantially less Hell absorption than the other two, as shown by the thin 
lines in Figure 3. We have rescaled the background intensity so that all three models produce the 
same mean absorption, but the greater emptiness of voids in the CCDM model remains evident in 
the density distribution of the absorbing material. 

When the same analysis is applied to HI absorption (thin lines in Figure 8), we see a shift 
towards higher density regions. The density regime pb/p~b < 2 accounts for ~ 55% of the HI flux 
decrement, compared to 80% for Hell. Regions below the mean density produce ~ 65% of the Hell 
absorption but only ~ 35% of the HI absorption. Small differences between the models appear 
mainly because scaling the UV background to match rneii does not give exactly the same thi in 
each case. 

The baryon density pb is not directly observable, but it is well correlated with the HI optical 
depth, which is observable. The upper left panel of Figure 9 is analogous to Figure 8, except that 
we use a threshold in thi instead of a threshold in pb, and we convert the mean flux decrement 
to an equivalent r in order to improve the dynamic range of the plot. We set the absorption 
of regions with thi below the threshold to zero before computing the mean decrement, D, and 
t = — log e (l — D). At high values of thi, the HI curves (thin lines) and Hell curves (thick lines) for 
each model converge, since saturated regions completely absorb both HI and Hell. However, much 
of the Hell absorption is associated with regions that have low HI optical depth. For example, half 
of the contribution to THeii comes from regions with thi < 0.15. Only ~ 15% of thi is produced 
by these regions. Figure 2 demonstrates that low density regions are revealed more clearly by Hell 
absorption than by HI absorption. Figure 9 demonstrates that these regions are in fact responsible 
for much of the Hell opacity of the high redshift universe. 

The upper right panel of Figure 9 compares the SCDM model to the two line models. The 
thermal and turbulent models are based on identical HI line populations, so their HI curves are 
identical, and their Hell curves are very similar. Even in the line models, much of the Hell 
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Fig. 9. — Effective mean optical depth as a function of thi threshold, at z = 2.33, analogous to 
Fig. 8. The value of r = — \og e (F) is calculated after setting F = 1 for points in the spectrum that 
have optical depth thi below the threshold value on the horizontal axis. Thick lines show THcii an d 
thin lines thi- Left hand panels compare the three CDM models; right hand panels compare the 
SCDM model to the line models with turbulent and thermal broadening. HI curves are identical 
for the two line models, since they assume identical HI line populations. Top panels show results 
for unsmoothed spectra. Middle panels show results after Hell and HI spectra are smoothed with 
a Gaussian filter of FWHM Av = 890 km s _1 , corresponding to 3 A in the Hell spectrum. Bottom 
panels show results when the Hell spectra are smoothed but the HI spectra are not. Only Hell 
curves appear in the bottom panels, since the HI curves would be identical to those in the top 
panels. 
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absorption occurs in regions of low HI optical depth — weak lines and the wings of strong lines. 
Relative to the CDM simulations, however, both line models produce a much larger fraction of 
Tjjeii in saturated regions. This difference corresponds to the visual impression one obtains by 
comparing sample spectra (Figures 4 and 7). 

The DKZ data have a spectral resolution of approximately 3A (Gaussian FWHM), so the 
predictions in the top panels of Figure 9 cannot be compared directly to the DKZ observations. 
The middle panels show the same analysis after the model HI and Hell spectra have been 
convolved with a Gaussian filter of 890 km s _1 FWHM, equivalent to 3A for Hell at z = 2.33. 
Regrettably, this smoothing eliminates most of the difference between the CDM models and the 
line models. The sense of this difference is reversed relative to the full resolution case, with the 
CDM models producing a slightly larger fraction of their absorption at high values of (smoothed) 
thi- This change probably reflects the large scale clustering that is present in the CDM models 
but not in the line models. Among the CDM models themselves, it is the more strongly clustered, 
CCDM model that produces the largest fraction of its absorption at high thi- 

The bottom panels of Figure 9 repeat the optical depth analysis using smoothed Hell spectra 
but unsmoothed HI spectra, in recognition of the fact that the best ground-based spectra are able 
to resolve even the narrowest observed HI Lya features. This version of the analysis preserves the 
clear distinction between the CDM models and the line models visible in the upper panels. Indeed, 
smoothing the line models amplifies the difference between them and the CDM models, since a 
significant fraction of the line models' smoothed Hell flux decrement arises in "inter-line" regions 
where the HI optical depth is extremely low. While a realistic comparison to observations will 
have to contend with noise in the Hell and (to a lesser extent) HI data, the distinction between 
the line models and the CDM models is likely to be observable. The distinctions among the CDM 
models themselves are more subtle. 

Figure 10 shows the distribution function of transmitted flux F; P(F)dF is the probability 
that a randomly selected point in the spectrum has transmitted flux in the infinitesimal range 
F — ► F + dF. Thick and thin lines represent Hell and HI results, respectively. Left hand panels 
compare the three CDM models, and right hand panels compare the SCDM model to the line 
models. Upper panels are based on unsmoothed spectra, lower panels on spectra smoothed over 
890 km s _1 . Figure 11 shows the distribution function of the log of the optical depth normalized 
so that J P(logr))dlog(r) = 1. While the information is equivalent in principle to that in 
Figure 10, the optical depth plot reveals tails of the distribution function more clearly. 

The distribution functions of the three CDM models are rather similar, though they are 
somewhat broader for the CCDM model because of its higher clustering amplitude. At full 
resolution, the Hell distribution functions of the CDM models are radically different from those 
of the line models, especially at low optical depths (F « 1). The line models have gaps in their 
spectra where the Hell absorption is extremely low, but the CDM models do not. The resulting 
difference in P(F) — an upturn at F « 1 for the line models but a downturn for the CDM 
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Fig. 10. — The distribution function of transmitted flux for z=2.33. P(F)dF is the probability 
that a randomly selected point in the spectrum has transmitted flux in the range F — ► F + dF. As 
in Fig. 9, thick lines show Hell and thin lines HI, left hand panels compare the CDM models, and 
right hand panels compare SCDM to the two line models. Top panels show results for unsmoothed 
spectra. Bottom panels show results after the spectra are smoothed with a Gaussian of FWHM 
Av = 890 km s _1 . 
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Fig. 11. — The distribution function of optical depth in logarithmic bins for z=2. 33. i- > (logr)(ilogr 
is the probability that a randomly selected point on the spectrum has the log of the optical depth 
in the range logr — ► logr + dlogr. The significance of panels and line types is the same as in 
Fig. 10. 
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models — also appears in the HI spectra, but it is less dramatic, and the downturn occurs at such 
low optical depth that it might be masked by errors in continuum fitting. The line models also 
have smaller fractions of the spectrum with transmission in the middle range 0.05 < F < 0.9. 
The -P(log r) distributions of the line models are much broader than those of the CDM models, 
reflecting the same trends seen in P(F). Unfortunately, the strong differences in the Hell spectra 
are mostly lost when the spectra are smoothed over 3A. The higher resolution spectra obtainable 
with HST may be better suited to detecting or ruling out the gaps in Hell absorption predicted 
by the line models, even though HST can only probe redshifts z > 3. Our simulations tend to 
underestimate the widths of P(F) and P(logr) for the smoothed spectra because the 890 km s _1 
filter width is a substantial fraction of our box size (2028 km s _1 at z = 2.33), and because we do 
not include spatial fluctuations in the intensity of the ionizing background (Zuo 1992; Fardal k, 
Shull 1993), which could be a significant source of additional fluctuations in the Hell opacity. A 
comparison of CDM simulations to ground-based, HI data using the cumulative form of P(F) is 
presented in Rauch et al. (1997) and Croft et al. (1997). For further discussion of this statistic, 
see Miralda-Escude et al. (1996). 

Figure 12 shows the distribution function of the log of the optical depth ratio, log(rH e n/THi)- 
As in Figure 9, the top panels show results for unsmoothed spectra, the middle panels show 
results with both spectra smoothed over 890 km s _1 , and the bottom panels show results using 
smoothed spectra for Tjjeii but unsmoothed spectra for thi- Starting with the top left panel, we 
see that the three CDM models predict similar distributions of log(rH e ii/THi), all peaked around 
THeii/iHi ~ 8. The extended tails towards low Tjjeii/ Thi are caused by thermal broadening, and the 
shorter tails towards high THeii/ Thi are caused by collisional ionization, as discussed in §3.3. In the 
smoothed spectra (middle left panel), the CCDM and OCDM models have broader distributions 
than the SCDM model, though the regime where this difference is strong is well below the peak 
of the distribution, and it may therefore be difficult to probe observationally. For smoothed Hell 
but unsmoothed HI (bottom left panel), the distributions become extremely broad and virtually 
identical. In particular, there are strong tails towards low THeii/THi caused by strong HI features 
whose Hell counterparts have been reduced by smoothing. The distributions also extend to high 
THeii/THi because the smoothed Hell optical depth is rarely less than 0.2, while the optical depth 
in the unsmoothed HI spectra can be 0.01 or smaller (see Figure 11). 

Right hand panels of Figure 12 compare the SCDM model to the line models. Without 
smoothing, the turbulent and thermal broadening models define two extremes, both very different 
from the CDM results. By construction, the turbulent broadening model has a constant ratio of 
THeii to thi, so its distribution is a 5-function at THeii/ Thi = 17. In the thermal broadening model, 
however, the HI and Hell 6-parameters differ by a factor of two, so in every line the ratio THeii/ Thi 
is high in the core and low in the wings. When the Hell and HI spectra are smoothed, the line 
model distributions are broader than the CDM distributions, especially in the direction of high 
THeii/THi- However, the strong differences again occur well below the peak of the distribution. 
Dividing smoothed Hell optical depths by unsmoothed HI optical depths again yields a broad 
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Fig. 12. — Distribution function of the log of the optical depth ratio, log(THeii/THi)> at z = 2.33. 
Layout is similar to that of Figure 9. 
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distribution function, but in this case the CDM results lie significantly above the line model results 
at the peak of the distribution, providing a clear distinction between the two sets of models. As in 
Figure 9, the difference reflects the presence of gaps in the line model HI spectra. 

Figures 9 and 12 show two useful 1-dimensional summaries of information contained in the 
2-dimensional, joint distribution of HI and Hell optical depths (Figure 5). Figure 11 shows the 
projections of this distribution along the HI and Hell axes. There are other potential cuts through 
this joint distribution, or one could use the full 2-dimensional distribution itself to distinguish 
between models. 

Our cosmological simulations have finite resolution, as discussed in §2. The gravitational 
softening length, 10 kpc comoving, is much smaller than the typical scale of Lya absorbers, but we 
do not adequately resolve baryonic structures less massive than ~ 32 SPH particles (4.65 x 1O 9 M 
for SCDM and CCDM, 2.11 x 1O 9 M for OCDM), and we cannot represent the initial fluctuation 
spectrum for wavelengths smaller than twice the initial particle grid spacing, A = 350 hr 1 kpc 
comoving. If we simulated the same volumes with much larger particle numbers, we would 
therefore expect some differences of detail for the statistical results illustrated in this section. 
However, the regions of the IGM that dominate Hell and HI absorption are low density and fairly 
smooth, and they lie outside the density/temperature regime in which cooling instabilities play an 
important role, so we do not expect major changes in our results to appear at higher resolution. 
We cannot firmly estimate the quantitative impact of resolution effects until we are able to perform 
simulations with more particles, which will be possible with a parallel TreeSPH code now under 
development (Dave, Dubinski & Hernquist 1997). We expect that the main qualitative effect of 
increased resolution on the simulated IGM will be a larger amount of low-amplitude substructure 
in underdense regions. Simulations with a factor of eight fewer particles than the ones used here 
yield similar physical properties for the IGM and a similar column density distribution for HI Lya 
lines with iVm <^ 10 15 cm -2 , but many fewer systems at higher column densities, where radiative 
cooling becomes important (Miralda-Escude et al., in preparation). 

The principal physical uncertainty in our simulations is the impact of reionization on gas 
temperatures. As pointed out by Miralda-Escude & Rees (1994), the energy injection of one 
photoelectron per proton during reionization can heat the IGM to several x 10 4 K, if reionization 
occurs rapidly enough that this energy is not dissipated by collisional line cooling. Our equilibrium 
treatment of photoionization suppresses this heating because we set neutral fractions to low 
values as soon as the ionizing background switches on, without altering gas temperatures. This 
treatment is equivalent to assuming that reionization occurred slowly (with consequent radiative 
cooling) or at high redshift (with consequent heat losses to Compton and adiabatic cooling). If 
we instead assumed that the IGM at z ~ 2 — 3.5 retained significant heat from reionization, then 
the temperatures of the unshocked or weakly shocked gas would be higher and less dependent on 
density. Recall that for the gas that produces most of the absorption, the HI and Hell optical 
depths are proportional to /?(,T _0 ' 7 r _1 . With hotter gas, the required values of Chi an d Ch c ii 
would therefore be higher, i.e., our models would require a less intense UV background or a 
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higher baryon density in order to match the observed thi and rneii- For a fixed reionization 
history, the (f^/i 2 ) 1 - 7 dependence in equation (8) would be closer to (f^/i 2 ) 2 , since the effect 
of reionization heating on gas temperatures does not depend on 1^ (though the subsequent 
photoionization heating does). The magnitude of these effects depends on the reionization model; 
plausible models could increase Chi and CHeii by 25 — 50% (see the discussion by Hui k, Gnedin 
1997). However, once the simulations were normalized to the observed tri and TR e n, we would 
expect the statistical properties of the absorption to be very similar to those computed here, with 
small changes reflecting the weaker dependence of optical depth on density and the larger degree 
of thermal broadening. Plausible changes in the reionization history all go in the direction of 
increasing rather than decreasing IGM temperatures, so they only tend to raise the lower bounds 
to flf, discussed in §3.2. IGM temperatures and corresponding flf, limits could be lower if the UV 
background spectrum is substantially softer than the HM spectrum, so that the mean residual 
photoelectron energy is lower. 



4. Conclusions 

Many previous discussions of Hell absorption have focused on distinguishing "line blanketing" 
in the Lya forest from "Gunn-Peterson" absorption by the IGM. Our cosmological simulations 
undermine the premise of this effort, for they suggest that Hell absorption and the low column 
density Lya forest both arise in diffuse, smoothly fluctuating, intergalactic gas. Local maxima in 
the optical depth can be identified as lines, but individual features do not, as a rule, correspond 
to compact structures that are sharply separated from their environment. While the Lya forest 
and the Hell flux decrement are both manifestations of the IGM, the high ratio of Hell ions to HI 
atoms does lead to an important systematic difference between helium and hydrogen absorption: 
Hell absorption is stronger in the mean, and much of it arises in underdense regions that have 
low HI optical depth. The general picture of the IGM presented here is similar to that in the 
other numerical simulation papers cited in the introduction, and it has much in common with the 
semi-analytic models developed by Bi (1993), Bi, Ge, & Fang (1995), Bi & Davidsen (1997), and 
Hui, Gnedin, & Zhang (1997). These semi-analytic models lead to similar qualitative conclusions 
about the properties of the gas producing Hell absorption (Davidsen, private communication; 
Davidsen et al., in preparation). 

Some of the more specific conclusions from our analysis are as follows: 

(1) The CDM models account for the observed relative values of tri (from PRS) and tr c r (from 
DKZ and HAR) if the UV background has the spectral shape predicted by HM. Large changes in 
the spectral shape (or in the observational estimates) would spoil this agreement. 

(2) These models account for the observed absolute values of tri and tr c r only if (a) the overall 
intensity of the background is lower than predicted by HM, or (b) the baryon density is higher 
than our assumed value of f^/i 2 = 0.0125. If we set the background intensity equal to the 
HM value, then the SCDM and OCDM models require Vt h h 2 « 0.023, and the CCDM model 



-32- 



requires VL^h 2 ss 0.038. We have analyzed one SCDM simulation with a higher baryon density, 
Ct^h 2 = 0.03125. If Tneii is held fixed, then this model produces stronger Hell absorption than 
the original SCDM model, as expected. Once both models are normalized to produce the same 
THeii, the statistical properties of their Hell absorption (e.g., the measures considered in §3.4) are 
virtually identical. 

(3) The CDM models naturally explain the observed evolution of rneii, in particular the factor of 
two drop in THeii between z ~ 3.3 (HAR) and z ~ 2.4 (DKZ), provided ThcII evolves at the rate 
calculated by HM. Change in the UV background does not play a major role in this evolution - 
Tiieii grows by only 20% between z = 3.3 and z = 2.4. The strong evolution of thcii is instead 
driven by the expansion of the universe, which lowers gas densities and spreads absorbing material 
over larger frequency ranges. The high Hell optical depth measured by HAR does not imply that 
Hell reionization occurred at z < 3.3. If the HAR upper limit of ThcII < 3.0 outside the Q0302-003 
ionization zone is correct, then the photoionizing background along this line of sight cannot be 
much softer than HM predict. 

(4) Most of the Hell opacity is produced by diffuse gas that follows a well defined relation between 
temperature and density. This relation has its origin in the competition between photoionization 
heating and adiabatic cooling. For gas that lies on this temperature-density relation, the Hell 
optical depth is a simple function of density, th c ii oc p^ 6 . To a first approximation, one can regard 
a Hell (or HI) absorption spectrum as a non-linear map of the gas density along the line of sight. 

(5) A significant fraction of the Hell absorption arises in regions with density contrast 5 < —0.5, 
for which linear perturbation theory will give inaccurate results. Analytic methods that treat this 
density regime more accurately, such as the Modified Zeldovich Approximation of Reisenegger 

& Miralda-Escude (1995), may provide useful guides to the behavior of Hell absorption in 
cosmological models, especially in light of point (4) above. 

(6) The three CDM models that we have investigated predict similar statistical properties of the 
Hell absorption (e.g., distribution functions of th c ii and thcIi/thi), once they are normalized to 
produce the same mean absorption. The CCDM model predicts somewhat broader distribution 
functions than the SCDM or OCDM models because of its higher mass fluctuation amplitude. 
Because this model has emptier voids, it also requires a lower UV background intensity and/or 
higher VL^h? to reproduce the observed Tneii, and it requires a slightly softer background spectrum 
to account simultaneously for thi and thcii- On the scales of our simulation, the three CDM 
models have power spectra of similar shapes, and their rms fluctuation amplitudes differ by less 
than a factor of two, so we do not yet know how the absorption results might change for much 
steeper or shallower power spectra or for models that have very different fluctuation amplitudes 
on these scales. 

(7) The CDM models have very different Hell absorption properties from our "line models," which 
assume that the Hell absorption is produced by randomly distributed, Voigt-profile lines with 
dN/dNm oc iVyj 1 ' 5 and a lower cutoff at Nm = 10 12 cm -2 . In particular, the line models have 
gaps in which the Hell absorption is very low, and a larger fraction of their Hell absorption arises 
in regions of high HI optical depth. The fluctuating IGM of the CDM models produces fluctuating 
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Hell absorption, but absorption-free regions (th c ii <S 0.05) are very rare. The line models require 
a softer UV background spectrum in order to produce the observed rneii , especially if the lines are 
thermally broadened. The statistical differences between the line models and the CDM models are 
greatly reduced when the spectra are smoothed, so for distinguishing these scenarios it is desirable 
to use Hell spectra with the highest resolution practical. If the minimum column density in the 
line models is pushed well below 10 12 cm' 2 , the differences from the CDM models are less striking, 
as weak lines overlap to produce a continuous, fluctuating background. 

(8) About half of the contribution to THeii comes from regions that have thi < 0.15. In the CDM 
models, most of the absorption in this regime is produced by a smoothly fluctuating IGM rather 
than the wings of strong absorption lines. The absorption is not uniform, but because the HI 
optical depth is low and the variations are gentle, much of the corresponding HI absorption could 
be inadvertently removed from optical quasar spectra in the process of continuum fitting. 

The CDM simulations provide a model of the high redshift IGM. Once it is normalized 
to the observed values of THeii and Tri, this model yields a number of testable predictions. 
First, the shape of the UV background spectrum should be close to that predicted by HM, with 
Thi/ThcH ss 100, for 2 <J z <J 3.3. It is difficult to measure the background shape precisely 
independent of an IGM model, but tests using metal line ratios (Songaila & Cowie 1995) or 
comparisons of the HI and Hell proximity effects (see Zheng Sz Davidsen 1995) might be able to 
identify strong departures from the HM spectral shape. If one adopts the intensity of the HM 
background as a lower limit, then our models also predict that the baryon density exceeds the 
"conventional" big bang nucleosynthesis estimate of O^/i 2 = 0.0125 by a factor ~ 1.5 — 3 (see Rauch 
et al. 1997 and Weinberg et al. 1997 for further discussion). The most direct and quantitative 
predictions of this IGM scenario are statistical properties of the absorption such as those shown 
in §3.4. In particular, observational analyses can test the predictions that a large fraction of the 
Hell opacity arises in regions of low HI optical depth and that there are few regions where THeii 
itself is very low. The DKZ spectrum offers the largest redshift range for statistical analyses of 
Hell absorption. However, the distinctive features of this IGM model are most stringently tested 
at high spectral resolution, so HST observations at z > 3 can also play an important role. An 
anecdotal but significant argument in favor of this IGM model is provided by HAR, who show that 
strong Hell absorption (th c ii > 1-3) persists in a region where a Keck HIRES spectrum reveals no 
HI Lya lines down to a detection threshold tri = 0.05. 

Careful comparison between the simulations and Hell observations should include realistic 
treatments of noise, instrumental resolution, the redshift range of the observations, and so forth. 
The three variants of the CDM model considered in this paper yield fairly similar predictions, and 
given the inevitable limitations of Hell measurements with existing instruments, we anticipate 
that HI observations will be a more powerful tool for distinguishing among them. However, Hell 
absorption provides a vital test of this general scenario for high redshift structure, checking one 
of its basic predictions in a regime that is almost inaccessible to other methods. In theories of 
structure formation like those studied here, underdense regions of the universe harbor a substantial 
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amount of diffuse baryonic material. If these theories are even roughly correct, then the recent 
studies of Hell absorption in quasar spectra have detected this material at high significance. More 
detailed analyses of existing absorption data can test this claim. Future instruments that can 
observe Hell absorption at high spectral resolution could yield a precise view of structure in the 
sparsest regions of the high-redshift universe. 
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